arXiv:1504.05949v3 [hep-ph] 21 Jul 2015 


Prepared for submission to JCAP 


The Electroweak Phase Transition in 
the Inert Doublet Model 


Nikita Blinov,Stefano Profumo,®’^ Tim Stefaniak®’^ 

“Department of Physics, University of California, Santa Cruz 
1156 High St, Santa Cruz, CA 95064 
^Santa Cruz Institute for Particle Physics, 

1156 High St, Santa Cruz, CA 95064 

“Theory Department, TRIUMF, 4004 Wesbrook Mall, Vancouver, BC V6T 2A3, Canada 
“^Department of Physics and Astronomy, University of British Columbia, Vancouver, 
BC V6T IZl, Canada 

E-mail: nblinov@triumf.ca, profumo@ucsc.edu, tistefan@ucsc.edu 

Abstract. We study the strength of a hrst-order electroweak phase transition in the 
Inert Doublet Model (IDM), where particle dark matter (DM) is comprised of the light¬ 
est neutral inert Higgs boson. We improve over previous studies in the description and 
treatment of the finite-temperature effective potential and of the electroweak phase tran¬ 
sition. We focus on a set of benchmark models inspired by the key mechanisms in the 
IDM leading to a viable dark matter particle candidate, and illustrate how to enhance 
the strength of the electroweak phase transition by adjusting the masses of the yet undis¬ 
covered IDM Higgs states. We argue that across a variety of DM masses, obtaining a 
strong enough first-order phase transition is a generic possibility in the IDM. We find 
that due to direct dark matter searches and collider constraints, a sufficiently strong 
transition and a thermal relic density matching the universal DM abundance is possible 
only in the Higgs funnel regime. 
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1 Introduction 

The simplest extension of the Standard Model (SM) that includes two SU{2) Higgs 
doublets is known as the inert Higgs Doublet Model (IDM). In the IDM, the extra doublet 
has no coupling to SM fermions and is odd under a postulated new Z 2 discrete symmetry, 
whereas all SM fields are Z 2 -even. Such symmetry makes the lightest Z 2 -odd particle 
(LOP) from the extra doublet stable and, thus, a potential weakly interacting massive 
particle (WIMP) dark matter candidate. The symmetry also eliminates numerous terms 
in the interaction Lagrangian of the model containing an odd number of extra “inert” 
scalars. 

The IDM was introduced originally as a possible generic scenario for electroweak 
symmetry breaking (EWSB) [1]. Only subsequently was it realized that the IDM natu¬ 
rally features a WIMP DM candidate lai, possibly providing a thermal relic density 
compatible with the inferred universal DM abundance. Numerous studies have subse¬ 
quently investigated the DM and collider phenomenology of the model (see, e.g., Ref. 11^ 

DEI). 

An additional early motivation to consider the IDM as an appealing augmenta¬ 
tion of the SM scalar structure was to allow for a relatively heavy SM-like Higgs while 
remaining compatible with constraints from electroweak precision observables, and with¬ 
out large fine tuning mm- Although this motivation has somewhat faded after the 
discovery of a SM-like Higgs boson at the LHC with a mass of ~ 125 GeV [III I12j . this 
important discovery decreases the number of free parameters in the theory by one, and 
places interesting and stringent constraints on the IDM phenomenology mm- 
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In the present study we are concerned with the nature of the electroweak phase 
transition (EWPT) in the IDM, and, specihcally, with determining which physical pa¬ 
rameters drive the strength of the phase transition, making it more or less strongly 
first-order, or second-order. This question is intimately related with the possibility to 
produce the observed baryon-antibaryon asymmetry in the Universe at the electroweak 
phase transition: A strongly first-order phase transition (in a quantitative sense we shall 
make clear below) is a necessary ingredient to (i) achieve the necessary out-of-equilibrium 
conditions, occurring on the boundary of broken and unbroken electroweak phase, and 
to (a) shield a baryon asymmetry captured in the broken electroweak phase region from 
sphaleron wash-out. 

While necessary, a strongly first-order phase transition is not a sufficient condition. 
The CP violating sources of the SM are known to be insufficient to generate the neces¬ 
sary asymmetry in the number density of baryons compared to antibaryons during the 
electroweak phase transition. The unbroken Z 2 symmetry in the IDM precludes any 
new source of CP violation, and thus this model per se cannot accommodate successful 
electroweak baryogenesis (EWBG). However, the IDM might be in effect a good ap¬ 
proximation at low energy of a broader construction that includes such additional CP 
violating sources at higher energies. Several suggestions of plausible effective higher¬ 
dimensional operators have been made in the literature MM- We will not discuss this 
aspect any more, as it falls outside the scope of this study. 

The nature and strength of the electroweak phase transition in the IDM has been 
subject of several studies, with increasingly refined treatment of the effective poten¬ 
tial HMD- For example, Ref. m utilized only the high-temperature form of the effec¬ 
tive potential without including the zero-temperature Coleman-Weinberg terms. These 
were then shown to be quantitatively important for the phase transition strength in 
Ref. [12], where the full one-loop effective potential was used. Alternative SU{2)i rep¬ 
resentations of the inert scalar were considered in Ref. |20| where it was argued that in 
general, higher representations are less successful in satisfying experimental and theo¬ 
retical constraints, thereby further motivating the study of the doublet case. 

With the exception of Ref. im, the primary focus has been on the Higgs funnel 
regime (described in more detail in Section]^. Indeed, we will confirm the findings 
of Refs. [181 US] that this is the only region of parameter space that can successfully 
saturate the DM abundance and provide a strong-enough hrst-order EWPT. In Ref. [2T| 
it was emphasized that the IDM can be useful for the EWPT even if the LOP provides 
only a sub-leading component of DM. 

In this work we go beyond previous studies by utilizing a state-of-the-art treatment 
of the finite-temperature and zero-temperature effective potential including renormal¬ 
ization group, daisy resummation improvements and one-loop model parameter deter¬ 
mination. As we discuss in great detail in what follows, strongly first-order EWPT in 
the IDM requires sizable quartic couplings that enhance quantum corrections to masses. 
This is important in the context of DM phenomenology since DM particle production in 
the early Universe often relies on resonance and threshold effects m- In addition, we 
also ensure that the phase transition completes by evaluating bubble nucleation rates. 
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Unlike previous studies which primarily utilized large numerical scans of the pa¬ 
rameter space, here we take an orthogonal approach: we restrict our attention to a 
few benchmark models, motivated by the requirement of having a viable dark matter 
particle candidate and representing different features in the DM phenomenology. Based 
on these benchmarks, we then discuss how the EWPT depends on the physical model 
parameters. We identify the key physical inputs that drive the phase transition to the 
interesting regime where it is strongly enough first-order to accommodate successful elec- 
troweak baryogenesis. We will see that in all but one case the demand for a strongly 
first-order EWPT is in tension with either the relic abundance requirement or with 
experimental probes. 

Our central hnding is that the main driver of the strength of the phase transition 
is the mass difference between the lightest inert scalar and the heavier scalars. Thus, we 
extend the results of Refs. [ITIIIH] to other regions of IDM parameter space. Por large 
enough mass splittings, but for light enough heavy scalars, we find a phase transition 
strength (as measured by the ratio Uc/To a-s we discuss in detail below) which increases 
with the mass splitting. 

The remainder of this paper is organized as follows: In Section we give a brief 
introduction to the IDM, thereby clarifying our conventions, discuss quantum and finite- 
temperature corrections to the effective potential, and outline the computation of the 
phase transition strength. The essential features of DM phenomenology are reviewed in 
SectionIn Section]^ we study the electroweak phase transition in several benchmark 
models motivated by the various DM scenarios available in the IDM. We conclude in 
Section [H 

2 Phase Transitions in the Inert Doublet Model (IDM) 

2.1 IDM at Tree-Level 

The IDM is a particular realization of the general type I Two Higgs Doublet Model 
(2HDM) (see, e.g.„ Ref. |22] for a review) which features an additional Z 2 symmetry. 
The SM doublet H is even under Z 2 , while the new (inert) doublet <I> is odd. If we take 
<I> to have hypercharge -1-1/2, the most general renormalizable potential consistent with 
these symmetries is then given by m- 


Uo = + A 2 |$|^ + A 3 |H| 2|$|2 + A4|Ht^>|2 + h + h.c. . 

( 2 . 1 ) 


Conventionally, within CP-conserving Higgs sectors, the physical states are decom¬ 
posed into CP-even and CP-odd scalars. One should keep in mind, however, that in 
the IDM there is no observable that can actually distinguish between the CP-even or 
CP-odd character of the inert Higgs bosons. In the absence of a vacuum expectation 
value (VEV) for <I>, the doublets decompose as 



( 2 . 2 ) 
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Below we consider the thermal evolution of the effective scalar potential in the early 
Universe. In general, spontaneous breaking of Z 2 can occur, in which case we must also 
include a VEV for the neutral component of which we will indicate with cj). 

The lightest Z 2 -odd particle is stable, and potentially provides a viable particle 
dark matter candidate. The Z 2 symmetry also forbids Yukawa couplings of <I> to SM 
fermions (assumed to be even under Z 2 ), which eliminates tree-level flavor-changing 
neutral currents. Either H 01 A can be the LOP, and since gauge interactions with 
SM states do not distinguish between the two, they are effectively equivalent from the 
standpoint of phenomenology. Below we will indicate the LOP as H, but all statements 
made with respect to DM phenomenology and the electroweak phase transition remain 
true after the replacements H ^ A and \l —>■ A^, where Xl = (A 3 -|- A 4 -|- X^)/2 and 
= (-^3 + A 4 — A 5)/2 determine the coupling of the LOP to the SM Higgs [T 8 ] . 

In the electroweak vacuum, the tree-level masses of the new states are given by 

m\ = n\ + 3Aiu^, 

m\ = /ii + ^ 

m-A = 

^‘h± = (2-3) 

The determination of model parameters from physical inputs is discussed in the next 
section. 


2.2 Finite-Temperature Corrections 

The effective potential at finite temperature T can be written as 

Heff = Uo + Hi + Vt, (2.4) 


where Vq, Vi and Vr are tree-level, one-loop temperature-independent and 
pieces, respectively. The tree-level potential Vq has been given in Eq, 


dependent 
(2.1). Working 


in the Landau gauge (.^ = 0), the temperature-independent one-loop correction has the 
Coleman-Weinberg form 


Hi = 


E Ui 

647r^ 


In 




-Q 


(2.5) 


The sum is over all particle species coupling to the doublets; rii is the number of de¬ 
grees of freedom (positive for bosons and negative for fermions), Ci are renormalization- 
scheme-dependent constants {Ci = 1/2 for transverse gauge bosons and 3/2 for every¬ 
thing else in the MS scheme); m/(u,(/>) is the field-dependent squared mass for each 
species. In writing the above, we have implicitly absorbed the counterterms into Ui; the 
temperature-dependent part is ultraviolet finite. The counterterms and, equivalently, 
the renormalized parameters and Ai are determined by ensuring that the one-loop 
potential reproduces the physical values v = 246.22 GeV and m/j ss 125 GeV. The re¬ 
maining parameters of the model are specified using the three physical masses niH, tua 
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and mH ±, along with A l and A 2 • The masses are related to potential parameters using 
the one-loop relations from Ref. m. while Xl and A 2 are taken to be running MS values 
defined at scale Mz- These requirements specify the renormalization conditions. 

The field dependent masses in the IDM for the SM vector bosons and fermions are, 
respectively, 

mw = ^9^iv^+ = + + = 0 ( 2 . 6 ) 


and 



- 2 2 


(2.7) 


with the corresponding bosonic degrees of freedom n, = 6, 3, 2 for i = W, Z, A, and 
fermionic degrees of freedom n, = —12, —12, —4 for i = t, b, r. 

The field-dependent neutral CP-even, CP-odd and charged scalar mass eigenstates 
are obtained by diagonalizing 


Ml = 
Ml = 
m| = 


+ Xl(I)^ 2Xl4>v \ 

(2.8) 

2Xict)v + ^X2ctP‘+ Xlv"^ J 

+ Am^ -b Xs^"^ X5cj)v \ 

(2.9) 

X5(j)V gL2 + + Xsv^) 

/if-b Ain^-b 5A3(/>2 i(A5-b A4)())n \ 

(2.10) 

l{X5 + X4)(l)v /if-b A2</>^-b jAsu^y ■ 


Equations (2.9) and (2.10) include contributions both from physical states and from 
Goldstone bosons. Notice that for 0 = 0 the (22) components reduce to the expressions 
in Eq. (2.3). 

The leading order quantum corrections give rise to a renormalization scale-dependent 
potential. One can choose the renormalization scale Q to minimize the size of higher 
order Ploop corrections which scale with (Inm^/Q^)^. The scale choice can be im¬ 
portant when a parameter in the potential is very different from the electroweak VEV 
~ 246 GeV. We thus choose to use the renormalization group (RG) improved effective 
potential to minimize the scale dependence. The potential parameters are replaced by 
their running values, evaluated at the scale Q. The relevant one-loop (3 functions are 
given in Appendix]^ Our computations are performed with Q = 246 GeV. 

The leading order temperature-dependent corrections to the effective potential in 
the Landau gauge take the form 


Vt = ^ ( X] ^ niJp [m^ (n, ()))/r^] j , 

\ 2 =bosons i=fermions / 

where the J functions are defined as 


Jb{x) = J dtt'^ In 1 — exp y—\/t‘^ + x 

poo p . 

Jf{x)= / dtt^ln 1 -|- exp (— + x 


( 2 . 11 ) 


( 2 . 12 ) 

(2.13) 
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These functions admit useful high-temperature expansions which allow us to study the 
phase structure as a function of T analytically (as long as the expansion is justified): 


m 


rpA J r 2 /rp2'\ rp2 2 / 2^3/2 1 4l 

r*[m/r]=-—+ -rm --T(m) 


Tvr^T^ 


T'^Jf = - - —T'^rn^ -m^ln ^ 

L / J 360 24 32 a/T2 


+ O (mVr^) 

(2.14) 

+ C>(mVr^), (2.15) 


where Of, = 16a/ = IGvr^ exp(3/2 — 27 ^). The terms in the expressions above illustrate 
symmetry restoration at high temperatures. The non-analytic term in Eq. (2.14) can 
be responsible for the barrier between the high T phase (at the field origin) and low T 
phase that breaks SU{2)l x U{1)y- 

Note that symmetry restoration signals the breakdown of perturbation theory — 
higher order diagrams become important. This can be accounted for by performing 
a resummation of daisy diagrams |26H28] . The resummation is performed by adding 


(2.16) 


finite-temperature corrections to the boson masses in Eq. (2.12): 

2 


m 


+ cT'^, 


where c is computed from the infrared limit of the corresponding two-point function. 
Eor the SM Higgs doublet we find 


Cl — 16^^^ 2”^^ 12'^'^ 12'^'^ 12^^ 4^* 4^^ 12^^’ 


The various components of the inert doublet receive similar contributions (but without 
contributions from the fermions): 


These expressions are in agreement with those in Refs. |191129( 130]. We implement these 


The thermal masses of the gauge bosons are more complicated. Only the longitu¬ 
dinal components receive corrections. The expressions for these in the SM can be found 
in Ref. |29j . but it is easy to modify them to include the contribution of an extra Higgs 
doublet. For the longitudinally polarized W boson, the result is 

~ + 2g^T‘^. (2-19) 


corrections by replacing /i? 

2 T 0 |p] 


fj,f + CiT^ in the scalar mass matrices, Eqs. (2.8 


2.9 


This includes contributions from gauge boson self-interactions, two Higgs doublets and 
all three fermion families. The masses of the longitudinal Z and A are determined by 
diagonalizing the matrix 


1 

4 


(u 2 + <^ 2 ) 


9 


-99 

99 ' 9 '^ 


+ 


f2g^T^ 0 \ 

V 0 2g'‘^T‘^) ■ 


( 2 . 20 ) 


^ There are subleading thermal corrections to off-diagonal self-energies suppressed by additional powers 
of coupling constants and VEVs which are usually neglected. 
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The eigenvalues can be written as 


where 


m 


Zl 


,A. = >| + (5^+5'^)r2±A 


1 


= ( -m| + (5^ + g'^)T^ - g^g‘^T‘^{v‘^ + (^^ + AT^). 


( 2 . 21 ) 

( 2 . 22 ) 


2.3 Electroweak Phase Transition (EWPT) 


Armed with the finite-temperature effective potential, we now proceed to study the 
structure of the EWPT. The key property we intend to investigate is the transition 
strength, which sets the baryon number wash-out rate inside a bubble of broken phase 
(for a recent review of electroweak baryogenesis, see, e.g., Ref. [M]). In order to suppress 
sphaleron wash-out in the regions of broken electroweak phase, the relevant condition is 
typically quantified by requiring that [32] 


Tc 


> 1 , 


(2.23) 


where Vc is the Higgs VEV at the critical temperature Tc, defined as the temperature at 
which the origin is degenerate with the electroweak-breaking vacuum. 

Note that it has been shown that this baryon number preservation condition (BNPC) 
is a quantity which is manifestly not gauge-invariant [33]. A gauge invariant BNPC can 
be however derived from the high-T expansion of the dimensionally reduced effective 
action and the critical temperature Tc must be obtained using the gauge invariant pre¬ 
scription of Ref. [33] , which employs expansions in powers of h of the potential and VEV. 
Near the critical temperature, 0{h) contributions to the potential are as important as 
the tree-level terms, so the h expansion fails. This is also why an all-orders ring diagram 
resummation discussed in Sec. 2.2 is needed. A consistent gauge-invariant method for 
implementing the ring resummation for the effective potential evaluated at the mini¬ 
mum was also demonstrated in Ref. [33]. We will be interested in studying tunneling 
and nucleation temperatures, which require the evaluation of the potential away from the 
minima. For this reason, below we employ the standard BNPC of Eq. (2.23) and use the 
full one-loop effective potential to study IDM phases. We will argue that our results do 
not depend strongly on the issues of gauge invariance. We leave the full gauge-invariant 
treatment of the IDM to future work. 

Finally let us note that the physical phase transition does not begin at Tc, but rather 
at a lower nucleation temperature T^, at which the bubble formation rate exceeds the 
Hubble expansion rate. If this rate is too slow, the false vacuum is metastable and the 
transition does not complete. We evaluate the nucleation temperature for a given model 
with a first-order phase transition using the CosmoTransitions package 


3 Dark Matter 

The requirement of a thermal relic abundance for the LOP matching the observed DM 
density in the Universe of = 0.1199±0.0022 [35], or at least of not over-producing 
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such density via thermal production (“subdominant IDM”, see, e.g., Ref. [2T]) naturally 
selects four distinct sectors of the model’s parameter space: 

1. a low mass regime, with a LOP mass, niH, well below half the observed SM-like 
Higgs mass, mn < m/j/2; 

2. a resonant or funnel region, mu ~ rnhj2, i.e., a mass range where LOP annihila¬ 
tion proceeds predominantly through quasi on-shell Higgs s-channel exchange; 

3. an intermediate mass regime, with a LOP mass of mhl2 <C mn ^ 500 GeV; 

4. a heavy mass regime, with a LOP mass between 500 GeV a few TeV. 

In the first case, the low mass regime, the DM pair-annihilation predominantly 
proceeds via the pair production of the heaviest kinematically accessible fermion (r 
leptons, b quarks) through h exchange. The lower the LOP mass, the larger the Xl^s 
couplings need to be in order to produce a large enough pair-annihilation cross section. 
The allowed mass values range down to values close to the classical Lee-Weinberg lower 
mass limit for WIMPs [36], for this class of models somewhere in between 3 and 4 GeV. 
Direct detection limits from XENONIO m probe such combinations of masses and 
couplings quite tightly, such that only a small mass window below 5 — 7 GeV remains]^ 

As the mass of the LOP approaches the resonant condition mn ~ mhf2, the res¬ 
onant Higgs exchange allows for much smaller values of the Xl^s couplings, and direct 
detection constraints can be readily evaded. The relevant mass window left uncon¬ 
strained by XENONIOO [39| and LUX [lOl has a width of approximately 10 — 15 GeV 
centered around mfi/2 |41j . 

The mass regions above and below the resonance mn = rnh/2 are actually slightly 
different from each other: Above the resonance, the pair production of WW* in the 
final state of DM pair annihilation processes becomes increasingly important, even if 
Xl,s = 0) because the four-point interaction through gauge couplings, independent of 
Xl,s, starts contributing significantly. As a result the values of Xl,s giving the “correct” 
relic density are pushed to increasingly (with LOP mass) large, negative values. 

For larger and larger LOP masses, the cross section for LOP pair annihilation to 
gauge bosons becomes very large, such that the thermal relic density is systematically 
below the universal dark matter density for any combination of model parameters. Bar¬ 
ring non-thermal production mechanisms, in this intermediate mass region the LOP 
cannot be the dominant dark matter constituent 

Finally, at about mn — 500 GeV, for Xl,s — 0 cancellations between scalar t- 
and u-channel exchange diagrams and the four-point interaction diagram alluded to 
above allow, again, for a sufficiently large thermal relic density. Such cancellations are 
suppressed by driving Xl,s away from zero. Thus, tuning Xl,s for increasing values 
of mu generally allows one to achieve the correct relic density for mass values from 
ruH ^ 500 GeV up into the multi-TeV range. This heavy mass regime of the IDM can 
be seen as a low energy effective theory of a composite dark sector |42j . 

■^The exact limit depends on the different possible choices of nucleon matrix elements, especially those 
connected with the strange quark content of nucleons | 38| . 
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Mh 

Ma 

Mh± 

Al 

A 2 

Tc 

Tn 

VciTc 

/^77 

BMl 

66 

300 

300 

1.07 X 10“^ 

0.01 

113.3 

110.3 

1.5 

0.90 

BM2 

200 

400 

400 

0.01 

0.01 

116.1 

113.7 

1.5 

0.93 

BM3 

5 

265 

265 

-6 X 10"3 

0.01 

118.2 

116.3 

1.3 

0.90 


Table 1. Input parameters for the three benchmark scenarios discussed in the text along with 
critical and nucleation temperatures, the transition strength and the signal strength for h —>■ 77 . 
The masses, given in GeV, are pole masses and the couplings are specified at Q = Mz- 
Temperatures are also given in GeV. 


4 Benchmark Models 

The benchmark models specified in Ref. m demonstrated various aspects of DM phe¬ 
nomenology and the possibility for the IDM to influence the h —?■ 77 rate. Unfortunately, 
none of the suggested scenarios exhibits a strongly first-order EW phase transition. 
In this Section, we identify alternate benchmark models which can potentially yield a 
strongly first-order EW phase transition, while having disparate properties for the light¬ 
est Z 2 -odd particle. All our benchmark models are compatible with constraints from 
Higgs collider bounds and rate measurements, which has been explicitly checked using 
the tools HiggsBounds [43l - li5] and HiggsSignals |l 6 ], where the model predictions have 
been calculated using a SARAH-generated SPheno version |47H50] . In the following dis¬ 
cussion we mostly focus on the interplay between the dark matter phenomenology and 
the strengths of the EWPT. 

Our key finding is that the requirement of a strongly first-order phase transition gen¬ 
erally leads to a large mass splitting between the LOP and the other scalars in the IDM. 
Our benchmark models are summarized in Tab. along with the corresponding criti¬ 
cal and nucleation temperatures, as well as phase transition strengths, as parametrized 
by the ratio Vc/T^- In each case the masses of the A and are chosen to ensure a 
strongly first-order phase transition. In Fig. we show the dependence of the transi¬ 
tion strength on these parameters. The lines corresponding to BMl and BM3 terminate 
where the potential develops a non-inert {4> ^ 0 ) vacuum hrst during thermal evolution. 
This vacuum can then either continuously evolve into the SM/inert ((/> = 0) vacuum at 
T = 0 or it can persist to low temperatures. In the latter case, the EWPT can occur 
in two steps. Such models are viable if the inert vacuum is deeper than the new one at 
T = 0 and both transitions complete {i.e., nucleation rate(s) are large enough). Two 
step electroweak phase transitions have been investigated in detail in Refs. [SUES]. In 
this work we consider only simple one step transitions, hence the truncation. Notice 
that in Ref. m the strength of the EW phase transition in models with multiple phase 
transition steps was always weaker, see Fig. 3 and 4 in Ref. |19j . 

First, let us consider model BMl, where LOP production in the early Universe is 
predominantly set by near-resonant s-channel Higgs exchange. This scenario has been 
recently examined in the context of phase transitions in Ref. m- Even more recently, 
it has also been suggested as a possible explanation m of the Fermi-LAT gamma- 
ray excess (see Ref. [M] and references therein). As discussed above, the on-resonance 
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requirement forces mn ~ but allows Xl to be small enough to be consistent with 

direct detection constraints. Here DM production does not rely on interactions with A 
or so their masses can be essentially chosen freely, as long as the resulting quartic 


couplings Aj (through Eq. (2.3)) satisfy perturbativity and constraints from electroweak 
precision observables (EWPO), which we check with 2HDMC [55]. In order to satisfy 


the BNPC of Eq. (2.23) one needs to increase the coupling of the new scalars to h, 
which, in turn increases the splitting of A and relative to H. For this and the 
following models we choose ruA = fo minimize the impact of splitting these states 
from H on the Peskin-Takeuchi T function |3| and to reduce the number of parameters. 
This assumption can be easily relaxed, but the results are qualitatively similar. This 
benchmark represents the only class of scenarios where the thermal LOP relic density 
(which we calculated with the micrOMEGAs code |56j l matches the observed DM universal 
density, and where the EW phase transition is strong-enough first order. 

When niA, ^ 340 GeV, the corresponding loop corrections to ruH are large 

and require < 0. This causes a second minimum to appear in the potential at T = 0. 
As ruA = is increased further, this minimum quickly becomes deeper than the SM 
one, corresponding to the termination of the blue curve in Fig. at ruA ~ 350 GeV. 
This behaviour was also observed in Ref. |19| . 

In this scenario, the LOP mass rriH has been chosen slightly above the kinematic 
threshold of the decay h ^ HH in order to evade constraints from direct searches for 
invisible Higgs decays and Higgs rate measurements. These however become important 
for our benchmark scenario BM3 (see below). 

The second benchmark BM2 in Tab. represents the intermediate mass regime. 
Here annihilation into gauge bosons is efficient and DM is generally under abundant, 
unless there is a cancellation among different amplitudes m- The cancellation depends, 
as indicated above, on how close Al ,5 0, i.e., on how degenerate the IDM Higgs sector 
is. In our benchmark, such a cancellation requires mn ~ ruA ~ with a maximum 

splitting of ~ 10 GeV. These small splittings lead to small couplings of the new states 
to h and therefore an insufficiently strong phase transition. Thus the phase transition 
requirement forces thermal relic DM to be under abundant. The observed DM density 
can be explained here, however, by invoking non-thermal production mechanisms {e.g., 
the decay of a heavy particle) or with the existence of additional DM particles {e.g., 
axions). The multitude of “non-standard” production mechanisms has been recently 
reviewed in Ref. m- 

The final benchmark model, BM3, belongs to the light-mass regime, and is another 
example that requires further ingredients to be fully consistent with the phenomenology 
of the DM sector. For < rnhl2, decays of the SM Higgs to invisible final states 
become possible, with a decay rate [3| 


T{h HH) = 


Svrm/i 


1 - 


2 “ 

mt 


1/2 


(4.1) 


Requiring consistency with the observed 95% G.L. upper limit on the branching frac¬ 
tion, BR(/i —)• HH) < 17% |58|, provides a strong constraint on the coupling Xi of 
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iriA = TaH± [GeV] 


Figure 1. Phase transition strength as a function of the heavier IDM scalar masses, taking 
rriA = 'rau±. The remaining parameters are chosen as in the benchmark models of Table 
which are shown by black dots. The lines for BMl and BM3 terminate where the inert doublet 
develops a non-zero vev, ^ 0, as described in the text. 


|Al| < 0.007, while a large value |Al| > 0.4 is required to sufficiently deplete the DM 
abundance |13j . These problems can be remedied by softly breaking the Z 2 , which would 
allow H to decay |59i I60j . As in the previous example, another explanation for DM is 
then needed. An alternative possibility is to provide the LOP with new annihilation 
modes, e.g., to new light vector bosons m, or a mechanism to dilute the thermal relic 
density, such as an episode of late entropy injection [62l [63] . 

DM phenomenology aside, it is again easy to get a strongly first-order phase tran¬ 
sition with a large mass splitting between H and A, We note that this scenario 

requires a significant tuning of parameters, because a small LOP mass requires near can¬ 
cellation of tree-level and loop contributions. For A^ > 0, this leads to negative values 
of which can result in the appearance of a new </) 7 ^ 0 minimum. 

In all three cases, the first-order transition is driven by the non-analytic 
terms (see Eq. 2.14) due to A and 77^, while the gauge boson contributions are not 
as important. This explains the common feature of large splittings between H and 
A, among the benchmark scenarios. These lead to large couplings between h and 
the new states, enhancing the size of thermal corrections. This appears to be a generic 
requirement for increasing the strength of the phase transition in the IDM. In particular, 
for the benchmark scenarios BMl, 2 and 3 we have (A 3 , A 4 ,A 5 ) « (3.3,—1.7,—1.5), 
(4.6,—2.3,—2.0) and (2.7,—1.4,— 1 . 2 ), respectively, at the scale Q = 246 GeV. These 
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coupling constants are large, but remain perturbative at energy scales of interest, which 
we checked using the RG equations in Appendix It is also important to emphasize 
that thermal corrections to the crucial terms from A and are not subject to 

gauge invariance issues that affect the gauge sector contributions. As a result, we expect 
these arguments to remain valid in the context of a fully gauge invariant treatment. This 
can be further tested in a toy model with all gauge coupling constants set to 0, thereby 
completely eliminating gauge dependence from the effective potentialj^ We checked that 
such a simplihed analysis gives quantitatively similar results for critical temperatures 
and transition strengths when the scalar couplings are large. 

The high-T expansion of the effective potential also provides a simple explanation 
for the shape of the curves in Fig.[^ In this limit the transition strength is proportional 
to the coefficient of the term [25]. For the IDM scalars such terms arise from the 
non-analytic contributions proportional to (assuming rriA = as 

above, and ignoring daisy contributions for simplicity), which behaves as only when 
;§> /i|. Thus when <C the transition strength is independent of IDM 

parameters, corresponding to the plateau of the green curve in Fig. In the opposite 
limit, the IDM gives an additional contribution to the cubic coefficient, so the transition 
strength scales as Vc/Tc ~ ^ m^, as illustrated by the monotonically increasing 

sections of the curves in Fig. |lp| 

For heavy masses w?/T'^ ^ 1 (with T ~ 100 GeV), the IDM states thermally 
decouple, but this does not mean that they have no impact on the phase transition. 
When /X 2 ;§> |/Uip, the heavy doublet can be integrated out to yield a SM effective 
theory with the potential 


Vo = + \\H\^ + + ... (4.2) 

where the dots stand for higher mass dimension operators. The parameters /r^, A and k 
can be related to those in the fundamental IDM by equating the effective potentials for 
the two models at a matching scale Q ~ /U 2 . For example, one-loop matching yields 

247r^^^" +^^ + ^ 3 / 4 ), (4.3) 

while fj? and A are determined below the matching scale by fixing the VEV and the 
Higgs mass. With the presence of a dimension-six term in the potential, the barrier 
required for a strongly first-order transition can be generated if A < 0 and -|- cT^ > 0 
for T ~ Tc, where c encodes thermal corrections from SM states only |29j . Scenarios of 
this type have been considered, e.g., in Refs. |64H67j . One immediate difficulty is that in 
the IDM K is generated only at one-loop, so in order for this operator to be significant 
for field values of around the electroweak VEV, one must overcome the loop suppression, 
suggesting that the combination A^-l-Al-l-Ag/d cannot be too small. This again forces a 

^We thank Michael Ramsey-Musolf for pointing this out to us. 

"^The precise scaling is modified by Daisy corrections, 0{n\/\svi) terms, £nite-T and renormalization 
group effects. 
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large splitting between the IDM states, meaning that the heavy DM scenario described 
in Section cannot be realized together with a strongly first-order phase transition. 
Such large couplings can run into perturbativity problems and invalidate the expansion 
used to generate the effective field theory. 

We briefly comment on the discovery prospects of the new IDM states at the LHC. 
Due to the Z 2 symmetry, the IDM states can only be produced pairwise at colliders. 
Successive decays of the heavier IDM states A and into the LOP and a Z or W boson, 
respectively, can give rise to multilepton signatures [TOl EH |69] . In a recent analysis [69] 
of LHC searches for supersymmetric particles with two leptons plus missing transverse 
energy in the final state in the context of the IDM, mass limits of up to ttia ^ 140 GeV 
for LOP masses mj/ < 55 GeV and charged Higgs masses around 85 — 150 GeV have 
been derived. While these limits partly exceed previous limits from the LEP collider, 
they are not yet sensitive to the parameter regions that yield a strongly hrst-order phase 
transition as required for successful electroweak baryogenesis, see Fig. 

The Higgs portal coupling can be probed directly at the LHC by searching 
for, e.g., jets plus missing energy m- However, the small magnitude of Xl in our 
benchmarks (as required by direct detection or relic abundance constraints) suggests 
that these searches will not be sensitive to the models in Tab. |TJ 

The new IDM states can also have an indirect effect on precision Higgs measure¬ 
ments. In particular, the new charged state provides an additional contribution to 
the loop-induced h —)• 77 and 7 Z rates. These effects have been recently studied in 
Refs. [711172j in the context of the 125 GeV Higgs boson. Modifications of these branch¬ 
ing fractions by 0 ( 10 %) are a generic feature of our benchmark scenarios, as we show 
below. The /i —)• 77 rate has the form |221|7IH75| 


r{h 77 ) 


a^GFm\ 

128v/27r3 


AIsm + 


Xzv"^ 

2 m|^± 


^0 




H± 


2 


(4.4) 


where the leading contributions to the SM amplitude AIsm ~ —6.56 -|- 0.08i come from 
W bosons and top quarks. The second term is the new contribution from 77^, where 
Ho is a loop function with the property lim 2 ;_s.o Ho(x) = 1/3 [75|- For our benchmarks 
we have A 3 > 0, so one expects a suppression of /i —>■ 77 relative to the SMj^Note that 
for fixed /r|, the amplitude for the contribution tends to a constant value 1/3 as 
A3 is increased. This means that in the limit of a large mass splitting between H and 
which is required for a strongly first order phase transition, the branching fraction is 
reduced by ~ 10%. This effect was also noticed in Refs. HEED for models similar to our 
BMl and BM2, respectively. The deviations to BR(/i —>■ 77) induced by are shown 
in Tab. [^in terms of the SM normalized signal strength = [a BR)/ {a BR)gM- While 
they are still consistent with the present measurements from ATLAS |76| and CMS [77| , 
the LHC should reach a precision of 4-8% for (58] EH] , thereby dehnitively testing 
the benchmark scenarios in Tab. [TJ 

®Various limits on LOP-Higgs coupling discussed above force |Ai| to be small, such that the H mass 
is primarily determined by fi 2 (at tree level, see Eq. ( |2.3[ )). If the charged Higgs is heavier than H 
then this forces A3 > 0. 
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While our benchmarks feature sizable deviations of BR(/i —)• 77 ) from the SM 
expectation, we note that it is possible to avoid this by taking to be nearly degenerate 
with H, and using A alone to drive the phase transition to be strongly first order. 
However, in this case, efficient coannihilation of H with during freeze-out generally 
results in a very small relic abundance m- The near degeneracy is also required by 
constraints on the oblique T parameter when niA ^ |3]- example, taking 

mjj± = 70 GeV, ruA = 370 GeV and other parameters as in BMl results in a strongly 
first order phase transition, an order of magnitude smaller relic abundance and only a 
~ 3% depression of relative to the SM. 

5 Discussion and Conclnsions 

We studied the structure of the electroweak phase transition in the inert Higgs doublet 
model, utilizing a set of three benchmark scenarios that feature a potentially viable dark 
matter particle. Our choices for the three benchmark models essentially exhaust all 
possible prototypical setups for particle dark matter in the inert doublet model. While 
only one of the benchmarks has a dark matter particle with a thermal relic density 
matching the observed dark matter density, the other two (under- and over-abundant) 
can be made viable by invoking additional production mechanisms or a scenario where 
the thermal relic density is diluted away, respectively. 

The key finding of our study is that in all cases where the model possesses a 
reasonable particle dark matter candidate, the inert scalar spectrum can be arranged in 
such a way so as to produce a strongly first-order electroweak phase transition. Central 
to achieving such a phase transition is to postulate a large enough splitting between the 
dark matter candidate and the heavier inert scalars. The physics driving this result is 
simple: Large mass splittings generically correspond to large couplings between the inert 
scalars and the Standard Model-like Higgs; These, in turn, increase the magnitude of 
non-analytic ~ terms in the temperature-dependent effective potential and thus 

the potential barrier between the field origin and the SU{2) x [/(l)-breaking phase. It is 
clear that this physical effect is generic, and, in fact, it has also been observed in other 
implementations of Two Higgs Doublet models |80H82j . 

The mass splitting under consideration cannot be arbitrarily large. For large enough 
values, for example, the phase structure of the model becomes more complicated, with 
possible non-zero vacuum expectation values for the inert doublet and multiple-step 
phase transitions. While, based upon the results of Ref. m the latter possibility is not 
expected to yield stronger electroweak phase transitions than in the single-step case, this 
is an interesting possibility which we leave for future studies. 

The question of how to embed large-enough CP violating sources in detail was 
also left unanswered here. It will be interesting to study whether such a source (for 
example an additional gauge-singlet complex scalar, see Ref. |83]) significantly impacts 
the electroweak phase transition and dark matter phenomenology, and, with this, the 
conclusions reached in the present study. 
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A Renormalization Gronp Equations 


Here we list the beta functions for the inert 2HDM at one-loop order. The general form 
of the RG equations is 


ca 

dt 


IGvr^ 


/3a, 


(A.l) 


where t = InQ/Qo and Qo is a reference scale. We take Qo = Mz- The U{1)y gauge 
coupling has the GUT normalization: gi = \Jb/2>g'. The beta functions below have been 
checked with SARAH m- Partial one loop results can be found in Refs. PI IM] for 
dimensionless parameters only. These agree with the formulae below. 

The gauge coupling evolution is determined by 


« 21 3 
^91 = y5i 

Pg2 = 

^93 =-7 al- 


(A.2) 

(A.3) 

(A.4) 


Eor the third generation Yukawas we have 

^7 2 ^2 o2 9330 2 /Ar\ 

Pyt = -^aiVt - |522/t - 8532 /i + + ^VbVt + VtVr (A.5) 

= -\alyb - ^ahb - ^alvb + ^ybVt + ybyl + ^jI ( a . 6 ) 

Pyr = -\alyT - ^ahr + "iahr + + ^yl- (A.7) 


Next we consider the scalar potential parameters. The evolution of the dimension¬ 
less quartic couplings A* is governed by 
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On n 27.900 

/^Ai = — — 9(72-^! + 

+ Qfl'2 + 24 A^ + 2 A| + A4 + A| + 2A3A4 

O 

+ I2\iyl - ^vt + l2Xiyl - Qyt + 4Aiy? - 2yt 

/3a 2 = -^51^2 - 95'|A2 + ^51 + + g52 + 24 A 2 

+ 2 A 3 + A 4 + Ag + 2 A 3 A 4 

9 o, ^ o, 27 4 9 o o 94 

/^Aa — — ^51^3 — 952^3 + “ ~^929i + ^52 

+ 4A3 + 2A4 + 2A5 + 12A1A3 + 12A2A3 + 4A1A4 

+ 4 A 2 A 4 + QX^yl + ^X^yl + 2X^y‘^ 

/3a4 = “^51^4 — 9(7|A4 + -9291 + 4 A 4 + 8 A| 

+ 4A1A4 + 4A2A4 + 8A3A4 + QXi^yf + 6X49^ + 2X49“^ 

Pxs = —-^ 9 i ^5 — ^ 92^5 + 4A1A5 + 4A2A5 + 8A3A5 

+ I2A4A5 + QX^y^ + GAsy^ + 2 X^y^ 

The beta functions for the mass parameters are given by 

- ^91^1 + 12Ai7r? + 4 A 3 //I + 2 X 49 I 
+ ^ylvt + ^ylvb + “^ylyr 

= ~Yq9i92 ~ 2^292 + 12A27r2 + 4 A 3//4 + ‘2‘X4y\. 

Finally, the anomalous dimensions for the Higgs and the inert scalar are 

9 2 92|o2,o2, 2 

Ih — ~'^9l — -^92 + 3% + 3yb + yr 

9 2 9 2 

' y < f > — “ 452- 


(A.8) 

(A.9) 

(A.IO) 

(A.ll) 

(A.12) 

(A.13) 

(A.14) 

(A.15) 

(A.16) 
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